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We apply Bogoliubov compensation principle to the problem of forward-backward asym- 
metry in tt production at TEVATRON. The non-trivial solution of compensation equa- 
tion for anomalous four-fermion interaction effective interaction of doublets of heavy and 
light quarks leads to possibility of existence of tu, bd vector bound states (resonance). 
The coupling of the bound state to t, u pair is calculated. With this result we obtain 
satisfactory description of totality of data on the effect of FB asymmetry including slope 
of effect dependence on invariant tt mass. Predictions for the resonance mass and pro- 
duction cross-section at LHC are presented.. 

Keywords: compensation equation; effective four-fermion interaction; forward-backward 
asymmetry. 

PACS numbers: 12.15.-y 11.30.Qc, fl.f5.Wx 

Recently results were presented on measurement of Forward-B ackwar d asym- 
metry in tt production in pp collisions at energy = 1960 GeV HEEl Afb = 
0.158 ± 0.075^, Afb = 0.196 ± 0.065 i^, ApB ^ 0.162 ± 0.047^1. According to 
the results the asymmetry significantly exceeds the corresponding SM calculations 
Afb = 0.089 ±0.007 SI. 

In present work we make attempt to explain the anomaly in terms of non- 
perturbative effects in SM, considering spontaneous generation of non-local four- 
quark interaction. In previous works we have used Bogoliubov compensation method 
as a mechanism for generation of non-local vertices, which produce non- 
perturbative effects in corresponding theories. Particularly, in work '-^ we obtained 
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Bogoliubov compensation equation solution for tt quarks' interaction, and made 
upon this ground some conclusions about composite Higgs properties. Below we 
shall start from this solution, but just for our present problem we have to take into 
consideration also mixed ut quarks' interaction. Accordingly, we have to investigate 
compensation equation for corresponding vertex and watch for influence, which it 
can produce in anomalous t-quark interactions. 

Bogoliubov compe nsation p rinciple (which quantum field theory adaptation 
in detail described in I9 | 10 | ll | allow us to obtain cutting form-factors for non- 
renormalizable vertices as a solution of special relations, expressing demand of 
mutual compensation of these vertices different order contributions into free quark 
Green functions. So, as a result we can use such free propagators in their usual form, 
and in the same time hold in our theory non-renormalizable interactions without 
divergence in integrals. 

Practically for mentioned relations construction we use "add-subtract" proce- 
dure: we insert additional term AL into Lagrangian both with "-|-" and "-" signs, 
and then one of this term we takes into free Lagrangian Lq, and other - into interac- 
tion Lagrangian Lint. So, we must then demand, that just the first addition together 
with constructed by this element higher order contributions into field propagation 
functions mutually compensate themselves. And the difference in sign in the same 
time will provide us preservation of analogous term in interaction Lagrangian: here 
compensation will not take place. 

In our theory we take following additional term: 



Here ^ and ■0 - correspondingly heavy quark doublet (t, b and light quark doublet 



with p,q,k,r - incoming or outgoing momenta of each line and F - form-factor func- 
tion. Firstly we consider special kinematic conditions, setting q = —p, k = r = 
as it presented in Fig.l and thus we have all diagrams depending only on one value 
p'^. Note, that color here is transmitted along lines. 



(1) 





P 



Fig. 1. Diagram for vertex ([T]). 
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The approximation, which we are working out here, proposes also following as- 
sumptions, where experience acquired in the course of previous works determines 
our choice: 

1) In compensation equation we restrict ourselves by terms with loop numbers 0, 1, 
2. 

2) We reduce thus obtained non-linear compensation equation to a linear integral 
equation. It means that in loop terms only one vertex contains the form-factor, 
being defined above, while other vertices are considered to be point-like. 

3) While evaluating diagrams with point-like vertices diverging integrals appear. 
Bearing in mind that as a result of the study we obtain form-factors decreasing 
at momentum infinity, we use an intermediate regularization by introducing UV 
cut-off A in the diverging integrals. Final results will not depend on the value of 
this cut-off. 

4) We shall not consider contribution to our equations of vertices of type 

^ ^^^L- Thus we multiply the equation by projecting operator 
(7p ~ just eliminating such terms and not affecting our standard vertices (see 
also El). 

These assumptions allow us to build compensation equation for vertex Gi , which 
in diagram form presented at Fig. 2 (as it was mentioned above, such equation for 
vertex G we obtained in previous work). 




+ 



= 



Fig. 2. Diagram representation of compensation equation for form- factor of 
interaction ([T|). Here solid lines represents t-quark, thin ones - u-quark. Also solid 
point represents vertex with form- factor F{p'^), whereas in the rest cases we have 
point-like four-fermion vertexes. 
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After angle integration we obtain following expression: 



3072 



^ ((10 G - 24 Gf ) In (x) x £f (y) dy 
24Gfln(a;)^ F {y) y dy + (^Hg, G-40Gf^ J F(y)ydy 



2G! + \ Gi g) i jy {y) y'dy + Q G, G - ^ Gf ) ^ ^ V (y) y^rf^^ 
24Gj2/ F (y) In (y) yd2; + :E (10 GiG- 24 Gf) / (y) In (y) dy 
(2 Gi G + 40 Gf ) a; / F (y) dy + (5 Gj G - 2 Gf ) / — ^dy 

^a; J X y 

^ G, G - ^ Gf ) £^dy + (36 Gf - 10 G, G) F (y) y dy 
(-10 Gj G + 24 Gf ) In (A^) +36 Gf ~ — Gj Gj x J F (y) dy 

/•OO /'OO \ 

+24 Gf In (A2) F (y) y dy + (30 G^ G - 48 Gf ) A^ F (y) dy j = . 

Here x = p^, y = with q - loop integrating momenta variable, and A - intervening 
momenta integration cut-off. Further results will not depend on this A. We count 
color number as iV = 3. 

Th is expre ssion provides an equation of the type which were studied in previous 
works l ^ i mril [ where the way of obtaining solutions of such equations was described. 
Indeed, by successive differentiation of this relation we obtain Meijer equation: 

(^4-) r.^ (^^ - 2^ r.^ - A (x± + 1] X 



dx J \ dx J \ dx J \ dx J \ dx J \ dx J \ dx 

^F(x)^ (144^- 12)a;2+ (^A.f (x)^ (8U - 48) x + F (x) (-30 C + 12)^ = 0. 
Here we take 

G = C G, . (4) 

Making substitution of the variable 

z = e — ; (5) 
we obtain solution in terms of a Meijer function 

F (2) = G GjJ (z Ij^g^ 1,1/2,1/2, 0; 0,~l/2,-l) ' (^) 



February 5, 2013 2:18 WSPC/INSTRUCTION FILE ttasymmetry 



Non-perturbative approach to anomalous FB asymmetry in top pair production 5 

where parameters ai, 02, 03, 04 depends on ^ and this choice of their distribution 
presuppose, that solution for ^ will give us two first, "active" parameters with 
negative real part, and two last - with positive. 

To determine ^ and C we consider two relation. First of them is normalization 
condition for form-factor F: 

^(0) = i; (7) 

Then we can mention, that we have inhomogeneous term with x — >■ of order of 
X In {x). So, it has to be eliminated by analogous term in Meijer function expansion. 
Near z = we have 



ffo2( \a^,a^;a,,a, \ , F (3/2 - Oi) F (3/2 Cz) jz j ( . .on 

^48 ^^13/2,1,1/2, 1/2,0; 0,-1/2,-1; T (03 - 1/2) F (04 - 1/2) " 



Taking also proportional x ln{x) term from the inhomogeneous part of the equation 
we obtain 

£^3/2-001^3/2-02) „ 
^'/' F(a3-l/2)F(o4-l/2) '^=V6 

Thus we obtain C, and substituting it into we have 

5 = 1.33623; (10) 

and, as a result, 

F (z) = 0.56229 G^? [z | 3;2"L'i72,i/2, o, -1/2, -1) 5 (H) 

where 

01 = -0.53192 -J 0.410037; 02 = -0.53192 + 1 0.410037 ; 

03 = 1.39470; 04 = 0.66913. (12) 

Since we have this solution for form-factor function, we can take into consideration 
bound state of interacting fields. Let us consider Bethe-Salpeter equation for a 
vector bound state. Its easy to see, that in our case this equation can be obtained 
by simple sign changing in front of the kernel of the compensation equation ([3]) 
(and also exclusion of inhomogeneous terms). So, it means, that we shall obtain as 
a solution Meijer function with the same parameters, as we just have in form-factor 
function, but with their changed distribution. Namely, we have such solution for 
normalized to 1 at zero momentum Bethe-Salpeter wave function ^ 

.T, I'-^N F(Q3)F(a4) ^42/ 01,02,03,04 \ 

* ^ r (1 - 01) F (1 - 02) V 3/2,1, 1/2,0,1/2,0,-1/2,-1^ ' U'JJ 



with ai (|12p . And then we can write normalization condition for bound state prop- 
agator, demanding, that linear by squared incoming momentum term in expansion 
of presented below diagram turn itself to 1. 
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t 




u 

Fig. 3. Diagram for normalization of the Bethe-Salpeter wave function. 



It gives value of coupling constant g of this bound state, standing in the vertices of 
diagram: 

9 = ^. h = ^ ; (14) 



/IT' ^ 2 

where we write fj, at lower integral limit instead of entering mass into propagator's 
denominator. We can estimate fi taking corresponding value fio from the results of 
heavy quark bound state investigation in work^^and accounting integration variable 
changing: 

/^o = 4.0675 * 10-l^ A^^^A^o (15) 

Thus we obtain following results: 

/i = 12.077, g== 0.904; (16) 

So, using these results we can estimate physical effects under discussion and 
compare obtained predictions with experimental data. 

Forward-backward asymmetry in proton collisions defined in such way: 

_ Njcosje) > 0) - N{cos{0) < 0) 

N{cos{e) > 0) + N{cos{0) < 0) ' ^ ' 

where 9 is the angle between proton and top-quark. CDF and DO Collaboration 
results for Aps are correspondingly 

Afb = 0.162 ±0.047 and Afb = 0.196 ± 0.065 . (18) 

We can also take into consideration invariant (t t) mass distribution of the asym- 
metry: 

Figure below shows CDF data together with SM predictions excluding background: 
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Parton-Level (GeV/c^) 

Fig. 4. Invariant mass M^^ distribution of the asymmetry. Lower Une presents SM 

calculations. 

At parton level we have the following SM prediction 

aM^j = (3.4 ± 1.2) X lG"^(GeV)"^ (20) 
and the experimental result ^ 

aM,- = (15.2 ± 5.0) X 10"'^(GeV)"i (21) 

Calculations of the effects are done with applications of CompHEP system 
We add to SM interactions our additional interaction between introduced bound 
state X and quarks: 

Lfx=9{^ 7a(l + 75)^ xX^+X^x + 75)*) . (22) 

For three values of parameter g from interval g — 0.9 ± 0.1 we can demand, that 
total cross-section of pp -> ti process reaches value of cr = 7.8 ± OApb and thus 
obtain corresponding value of M = Mx- Then accounting all possible sub-processes 
we can estimate for each g value of asymmetry A as also slope a. Results presented 
in Tables 1 and 2. 

Table 1. 

FB asymmetry for different values of coupling g and bound state X mass M 

(GeV). 



(7 


8.2pb 


7.8pb 


7Apb 


g 


M 


Afb% 


M 


Afb% 


M 


Afb% 


0.8 


528 


22.21 ±0.07 


550 


18.4 ±0.05 


575 


14.5 ±0.04 


0.9 


608 


21.4 ±0.06 


632 


17.8 ±0.05 


660 


14.1 ±0.04 


1.0 


686 


20.7 ±0.06 


713 


17.3 ±0.04 


747 


13.5 ±0.04 
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Table 2. 

Slope a in 10^''/GeV for different values of coupling g and bound state X mass M 

in GeV. 





8.2pb 


7.8p& 


l.Aph 


g 


M 


a 


M 


a 


M 


a 


0.8 


528 


18.23 


550 


17.14 


575 


16.34 


0.9 


608 


19.91 


632 


18.62 


660 


16.99 


1.0 


686 


19.41 


713 


17.33 


747 


16.32 



The results, presented in Tables, show, that our proposal can describe the ex- 
isting data on the effect. Indeed the calculated value s of asymmetry A and of slope 
a are inside error bars of experimental results l ^ l ^ l '^ l for whole intervals of variation 
0.8 < 5 < 1.0 and l.Apb <a < 8.2pb. 

It is interesting to compare with data from LHC. First of all the forward- 
backward asymmetry can not be measured at LHC due to symmetry of the ini- 
tial state pp. We also are to emphasize, that the dominant contribution to the tt 
production at LHC is given by gluon-gluon fusion, in which contribution of the 
mechanism under discussion is absent. However the production of the new states 
X, X can contribute to the total cross-section of tt production. The estimate of this 
contribution for g = 0.9 and M = 632 GeV gives the following value for ^/s = 7TeV 

Act = ll.Gpb; (23) 

where we have taken into account value of the branching ratio 

BR{X ->tu) ^ 0.A7; Tx = IMGeV . (24) 

Value (125)) is inside error bar s of t-p air production cross-section at LHC |14 | 15 | 16 | 
For example, results of works! I I read correspondingly 

cr{pp^ti+ ...) ^ 176±5tlt±8pb; (25) 
(T{pp~>tt + ...) = 158±2.1± 10.2±3.5p6; 

while the SM calculations give 

a{pp^ti+ ...) ^ 165tllpb. (26) 

We see, that possible additional contribution ()23|) causes no trouble yet. 

However we are to emphasize, that our variant predicts the production of X, X 
resonances and it is quite advisable to look for this effect in the forthcoming LHC 
studies. Let us remind parameters of these resonances: mass in the range 600 — 
700 GeV, width around 150 GeV and production cross section for y/s = 7TeV 
in channel X tu is oi the order of magnitude (|25)) . The resonance can be also 
observed in the channel X bd with branching ratio 0.53. Note, that our approach 
also give results for processes with b b production, which are of the same order of 
magnitude as for t-quarks. However, it seems, that background conditions for b b 
processes are more severe. 
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In case of confirmation of the present option we would get the result, that the 
asymmetry under discussion by no means contradict the Standard Model. On the 
other hand this would mean necessity to take into account non-perturbative effects 
in SM, in particular, by application of the Bogoliubov compensation approach. 
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